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We consider the topological susceptibility for an SU(N) gauge theory in the limit 
of a large number of colors, N —> oo. At nonzero temperature, the behavior of 
the topological susceptibility depends upon the order of the deconfining phase 
transition. The most interesting possibility is if the deconfining transition, at 
T = Tjj, is of second order. Then we argue that Witten’s relation implies that 
the topological susceptibility vanishes in a calculable fashion at Tj. As noted by 
Witten, this implies that for sufficiently light quark masses, metastable states which 
act like regions of nonzero d — parity odd bubbles — can arise at temperatures 
just below T,i. Experimentally, parity odd bubbles have dramatic signatures: the 
g' meson, and especially the g meson, become light, and are copiously produced. 
Further, in parity odd bubbles, processes which are normally forbidden, such as 
g —> TT^n^, are allowed. The most direct way to detect parity violation is by 
measuring a parity odd global asymmetry for charged pions, which we define. 


1 Introduction 

In this paper we give a pedagogical introduction to recent work of oursS We 
consider an SU{N) gauge theory in the lirwt of a large number of colors, 
N oo. This is, of course, a familiar limitEl We use the large N expansion 
to investigate the behavior of the theory at nonzero temperature, especially 
for the topological susceptibility. The results depend crucially upon the order 
of the deconfining phase transition; if it is first order, nothing very interest¬ 
ing happens. If the deconfining transition is of second order, however, the 
topological susceptibility vanishes in a calculable fashion. 

This implies that metastable states, which act like regions with nonzero 
0, can appear. Parity is spontaneously broken in such parity odd bubbles, 
and produces novel physics. The rj meson becomes very light, at most a few 
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hundred MeV, and so is easily produced. As parity is broken, the rj can decay 
into two pions, instead of the usual three. We also propose a global variable 
which can be used to measure an asymmetry in parity. 

2 Large N 

Holding the number of fermion flavors fixed as N oo, the large N limit 
is very much a gluonic limit, as the ^ iV^ gluons totally dominate the ^ N 
quarks. This is the basis for most of the conclusions which we can draw at 
large N: given what happens to the gluons, what happens to the quarks follows 
almost immediately. 

The standard assumptions at large N is that the physics is like that for 
N = 3: confinement occurs, so at low temperature we can speak entirely 
of mesons and glueballs. Their masses are assumed to be of order one as 
N —>■ oo] as usually occurs in any large N limit, interactions between either 
mesons and/or glueballs are suppressed by powers of 1/N. 

It is also natural to assume that the degeneracy of mesons and glueballs is 
of order one. This is, after all, what we mean by confinement: all trace of the 
color indices disappear, leaving bound states which are characterized only by 
spin, parity, etc. Thus in the low temperature phase at a temperature T, since 
each meson or glueball has a free energy which is ~ the total free energy 

is also of order one. In contrast, in the deconfined phase atiiigh temperature, 
the free energy is of order ~ N'^. As pointed out by Thornl^ this allows us to 
use the free energy itself as an order parameter for the phase transition: the 
transitiop._at a temperature T^, occurs when the term in the free energy ~ N'^ 
turns onffil 

Rigorously, the true order parameter for the deconfining phase transition is 
associated with the spontaneous breaking of the global Z{N) symmetry above 
Td] this symmetry becomes 0(2) as N ^ oo. 

We also make a further assumption, namely that any other phase transi¬ 
tions occur at the same time as deconfinement, at T^. We have no proof of this 
statement, although we suspect that a proof can be constructed in the limit 
of large N. Nevertheless, it strains credulity to image that at the point where 
this huge increase in the free energy occurs, that that alone doesn’t force any 
other phase transitions in the theory. 

In particular, assume that we couple massless quarks to the gluons. At 
zero temperature, it is known that the quarks’ chiral syrnmetry must be broken 
in the familiar pattern, to a diagonal subgroup of flavorH Then we assume that 
the chiral symmetry is restored at a temperature with = Td- 

We take the scale of the deconfining transition to be the same as for the 
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glueball masses; thus Td is of order one as iV ^ oo. This turns out to be 
a remarkably powerful assumption. Consider the large N limit of a theory 
without confinement, such as a A^-component vector with coupling 5 ^, holding 
g^N fixed as N ^ 00 . We assume that the masses of the fields are of order 
one. Then the only way for a transition to occur in what is, after all, free field 
theory, is to go to temperatures which grow with N ; a simple one loop estimate 
gives ~ ^1 ~ (this is also the scale of /^, which is natural). What 

happens in a confining theory is far more dramatic: the transition occurs at 
temperatures of order one, not ~ . This implies that the hadronic phase 

is “cold” at large iV: interactions are small, so that effects from the thermal 
bath, such as the loss of manifest Lorentz invariance, can be neglected. 

The crucial thing which we do not know about the large limit is the order 
of the deconfining phase transition. (The effect of quarks can be neglected, 
since the gluons dominate the free energy above Td\ this will be elaborated 
later.) For this we must look to the lattice, which as always provides the true 
intellectual basis for our understanding. 

In the early days of Monte Carlo simulations on the latticeP it was generally 
agreed that the deconfining phase transition is of first order when = 4. It is 
not clear, however, if these simulations are definitive. In particular, they were 
done at rit = 4, where rit is the number of steps in the imaginary time direction. 
For the standard Wilson action, at this value of rit there is a bulk transition 
close to the finite temperature transition. The bulk transition can be avoided 
by going to larger values of rij, but at the time this was computationally difficult 
to do. Recently, however, Ohta and Wingatal have computed for A^ = 4 and 
Tit = 6 ; they find that the strong first order transition at n* = 4 is gone for 
rit = 6 . Of course, to really establish that there is a true second order phase 
transition is a difficult matter, requiring lengthy study. But these results do 
suggest that it may be hasty to conclude from rit = 4 that the deconfining 
phase transition is of first order. 

There are also results on the large N limit of gauge theories on the lattice! 
Such reduced models appear to reliably predict the ratio of the critical tem¬ 
perature to the square root of the string tension. They predict a first order 
transition, but only under the technical assumption that the coupling between 
spacelike plaquettes can be neglected. It is not apparent to us how strong this 
assumption is. 

Previously, Pisarski and Tytga0 suggested that the large N deconfining 
phase transition is of second order. Their argument was rudimentary: the 
easiest way to understand why the deconfining transition is weakly first order 
for A^ = 3 is if the large N expansion is a good approximation, and if the 
transition is of second order for N = cx). Then the cubic invariant, which 
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drives the transition first order at iV = 3, is suppressed by ^ 1/iV. Of course 
the first assumption is rather strong: perhaps the large N expansion is not a 
good guide to thermodynamic properties. 

In the following we assume that the deconfining transition is of second 
order for all iV > 4, but comment upon how our results change if the transition 
is of first order. irn 

The principal object we are interested in is the topological suceptibilityU’Eil 
From the topological charge density, 


Qix) = {g^/327r^)tr{G^0G^^) = 


( 1 ) 


The current is gauge dependent. The topological susceptibility is the two 
point function of Q, 


Avm(T) = d^Fi9,T)/de^ = J d^xQ{x)QiO) ; 


( 2 ) 


F{9,T) is the free energy, and the 9 parameter is conjugate to Q. At zero 
temperature, the free energy reduces to the energy, F{9, 0) = E{9). 

Since Q is a total derivative, Xym{T) vanishes order by order in pertur¬ 
bation theory. It receives contributions entirely from nonperturbative effects, 
such as instantons. The action of a single instanton with fixed scale size is 
In the large N limit, g'^N is held fixed as N ^ oo, so the contribution 
of an instanton to the topological susceptibility is \ym{T) ~ exp{—aN), with 
a = Stt'^ / {g^^ N). Thus the contribution of instantons vanishes exponentially 
in the large N limit. This naive argument assumes that the integral over in¬ 
stanton scale size is well behaved. This is certainly true in the limit of high 
temperature; then the theory is weakly coupled, and instantons are suppressed 
by the Debye screening of electric fluctuations. This naive picture was veri¬ 
fied, at all temperatures, by Affle^ in a soluble asymptotically free theory, the 
GP^ model in 1 -|- 1 dimensionfiJ. 

Thus at large N, in the deconfined phase the topological susceptibility is 
exponentially small in 1/A^, and so essentially vanishes. At zero temperature, 
Witteiidl suggested that instead of semiclassical fluctuations, that quantum 
fluctuations generate a nonzero value for AyM(O) ^ A^°Edl It is natural to 
assume that the topological susceptibility is ^ N° throughout the deconfined 
phase, and changes to ^ exp{—o^) at the deconfining phase transition. This 
was previously argued by AffleclcJ and by Davis and MathesonEn 

How it changes depends upon the order of the phase transition. If the 
deconfining transition is of first order, then as the hadronic phase is “cold”, 
the most natural possibility is that the topological susceptibility is essentially 
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constant in the hadmnic phase, and changes discontinuously to zero at T^. 
Recent lattice resulted in full QCD with four flavors indicate a sharp drop for 
the topologocal susceptibility across the phase transition, and thus seem to 
support this conjecture. 

If the d£Confining transition is of second order, a more extended analysis 
is necessaryld Generalizing the results of WitteiEdI and Venezianetd to ncpzero 
temperature, and using results on the anomalous couplings of mesonsl^ we 
find that the free energy depends upon 9 as 

(l + c02)(T,-rf-“ . (3) 

Here a is the critical exponent for the deconfining phase transition, a « —.013. 
We then use Witten’sdJ formula for the 77 ' mass to conclude that the rj' mass 
vanishes at Td, 

= {Td-Tf-^. (4) 

Implicitly, we have used the fact that the hadronic phase is cold, so that zero 
temperature formulas, such as (i, generalize trivially. 

3 Parity odd bubbles 

We now use this result on the r]' mass to investigate the nature of the theory 
in the hadronic phase, just below Td- At zero temperature, a successful phe-. 
nomenology of the rj' was developed with a chiral lagrangian formalismEdllld 
For Nf flavors, a U{Nf) matrix U is introduced, satisfying U^U = 1. U de¬ 
scribes the Nj — 1 pions and the 77 '. The effects of the anomaly are represented 
solely by a mass term for the 77 ', {tr In U)^. 

This is stark contrast to how effects of the anomaly due to instantons are 
included. Consider a linear sigma model with a field $. Then the effects of the 
anomaly enter exclusively through a term ^ det{^). As in the nonlinear sigma 
model, when the chiral symmetry is spontaneously broken, this term generates 
a mass for the rj'. However, in the nonlinear sigma model, at large N, there is 
only a mass term for the 77 '; four point interactions between 77 '’s are induced 
by the anOjCaaly, but are suppressed by higher powers of 1/N'^. As emphasized 
by Witten^ a term ^ det{^) violates this large N counting. This is subject 
to the trivial qualification that Nf > 4, so that there are quartic interactions 
between the 77 '’s. 

For U fields which are constant in spacetime, the potential for U is 

f2 

V{U) = ^ {tr {M{U + C/t)) - a{tr In U - Of) ; (5) 
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The pion decay constant /tt = 93 MeV, while M is the quark mass matrix. 
When M = 0, ^ a, so a ^ /N. 

Taking Mij = , any vacuum expectation value (v.e.v) of U can be 

assumed to be diagonal, Uij = . The potential reduces to 

This is minimized for 

sin{(j)i) + a(Y^ (j), - 9) = 0 . (7) 

Note that as ^ 4>i arises from tr In U, it is defined modulo 27r. 

Previous^jSeveral authors studied how the v.e.v.’s of the (j)’s change as a 
function of SEjiU In the present work, we confer 0 = 0, but consider how 
the fXi and a change with temperature. WitteitJ pointed out that when the 
anomaly term a becomes small, metastable states in the (jj’s can arise. From 
our arguments in the previous section, this happens naturally if the phase 
transition at large N is of second order. 

The presence of these metastahle states can be easily undepikyd for a 
single flavor, as discussed by WitteicJ (for a recent discussion sefi^O). From 
(0). the v.e.v. arises from a balance between a term ~ sin{4>) and a term 
~ (j). For large a, the term linear in (j) wins, and there is no possibility for a 
metastable point. Now consider the opposite limit, of vanishing a: then there 
automatically other solutions besides (j) = 0, (p = 2 tt, Att, etc. These solutions 
are equivalent to the trivial vacuum, and so there is nothing new. But for 
small values of a, the term linear in a will only move the stationary point a 
little bit from 27r, Att, etc. Because a is nonzero, they will become metastable, 
distinct from the usual vacuum. 

From these states will act like regions of nonzero 9. Parity and CP 
are both violated spontaneously in such a region. 

The condition for metastable states to arisOj-ivith more than one flavor is 
not apparent, and a new result of our analysisEJ It is easiest understood by 
analogy. At zero temperature, and nonzero 9, if any quark mass vanishes, the 
9 parameter can be eliminated by a chiral rotation through that quark flavor. 
Thus it is the lightest quark mass which controls 9 dependence. We found a 
similar phenomenon for metastable states: they only occur when the anomaly 
term is small relative to the lightest quark masses. 

This means that metastable states only arise when the anomaly term be¬ 
comes very small. At zero temperature, the anomaly term is on the order of 
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the strange quark mass. The previous argument indicates that it must become 
on the order of the up and down quark masses. Putting in the numbers, we 
find that metastable states only arise when the anomaly term becomes on the 
order of 1% of its value at zero temperature. Clearly this is a strong variation 
of the topological susceptibility with temperature; nevertheless, it is interesting 
to investigate the possible implications for phenomenology. 

Most notably, when— the^nomaly term a becomes small, there is maxi¬ 
mal violation of isospinEjO’OO At zero temperature, the nonet of pseudo- 
Goldstone bosons — the it’s, AT’s, rj, and 77', are, to a good approximation, 
eigenstates of SU{3) flavor. It is not often appreciated, but this is really due 
to the fact that the anomaly term is large, splitting off the rj' to be entirely 
an SU{3) singlet. When the anomaly term becomes small, however, while the 
charged pseudo-Goldstone bosons remain approximate eigenstates of flavor, 
the neutral ones do not. Without the anomaly, the 7r° becomes pure uu, the rj 
pure dd, and the rj' pure ss. Consequently, these three mesons become light. 
This is especially pronounced for the 77, as it sheds all of its strangeness, to 
become purely dd. Thus the 77 and 77' would be produced copiously, and would 
manifest itself in at least two ways. First, light 77’s and rj"s decay into two 
photons, and so produce an excess at low momentum. Secondly, tb£se mesons 
decay into pions, which would be seen in Bose-Einstein correlation^^. Further, 
through Dalitz decays, the entp^ced production of 77’s and 77'’s will enhance 
the yield of low mass dileptons E3. 

This maximal violation of isospin is true whenever the anomaly term be¬ 
comes small. There are other signals which only appear when parity odd 
bubbles are produced. Since parity is spontaneously violated in such a bub¬ 
ble, various decays, not allowed in the parity symmetric vacuum, are possible. 
Most notably, the 77 can decay not just to three pions, as at zero temperature, 
but to two pions. Because of the kinematics, in a parity odd bubble, 77 ^ 
is allowed, but rj tt + tt " is not. 

There is another measure of how parity may be violated. We first argue 
by analogy. Gonsider propagation in a background magnetic field. As charged 
particles propagate in the magnetic field, those with positive charge are bent 
one way, and those with negative charge, the other. This could be observed by 
measuring the following variable globally, on an event-by-event basis: 


^=E 

TT+TT” 


fe+ X Ptt-) • ^ , 

\PTr+ I IP-k- I 


( 8 ) 


here z is the beam axis, and p are the three momenta of the pions. 

If the quarks were propagating through a backgmund chromo-magnetic 
field, then V, which is like handedness in jet physicsl^ is precisely the right 


7 




quantity. However, a parity odd bubble is not directly analogous to a back¬ 
ground chromo-magnetic field: tt+’s and 7r“’s propagate in a region with con¬ 
stant but nonzero (j) in the same fashion. Consider, however, the edge of the 
parity odd bubble: in such a region, Wdfj,U is nonzero, and does rotate 
and TT~ in opposite directions. Thus it is the edges of parity odd bubbles which 
contribute to the parity odd asymmetry of (H). Purely on geometric grounds, 
this suggests that a reasonable estimate for the maximal value of V is on the 
order of a few percent. 

We conclude by noting that what appears to be a rather technical subject 
— the 0 dependence of the free energy — is related to interesting and novel ex¬ 
perimental signatures in heavy ion collisions. Within our assumptions, we find 
that parity odd bubbles only arise very near the point of the phase transition. 
This is very much tied to the fact that we limit ourselves to an analysis at large 
N. For finite N, it is a long standing question of how to reconcile the known 
limit at large N with periodicity in 6, with period 2tt. A jprpbable solution 
involves “glued” potentials, which are a sum of cosines; se^’tZl and recentlyfis 
The precise form of the potential at finite N could jiramatically alter our re¬ 
sults, and, as discussed by Halperin and Zhitnitsky^ make the emergence of 
parity odd bubbles far more likely. 
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